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We discuss a class of phantom (p < —g) cosmological models. Except for phantom we admit 
various forms of standard types of matter and discuss the problem of singularities for these cosmolo- 
gies. The singularities are different from these of standard matter cosmology since they appear for 
infinite values of the scale factor. We also find an interesting relation between the phantom models 
and standard matter models which is like the duality symmetry of string cosmology. 
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I. INTRODUCTION 

The observations of distant supernovae has given 
the evidence for accelerating universe filled with matter 
which violates the strong energy condition g + 3p > 0, 
where p is the pressure and g is the energy density. 
However, more detailed check of data suggests that also 
the matter which violates the weak energy condition 
p > OjP > is admissible at high confidence level 

The physical background for strongly negative pres- 
sure matter (phantom) may be looked for in string the- 
ory . The point is that the group velocity of a wave 
packet is negative which leads to the behaviour of the 
packet which is somewhat unusual - namely, it moves in 
the direction opposite to the momentum. In such circum- 
stances a wave packet which leaves the comoving volume 
transfers momentum into the volume which makes nega- 
tive contribution to the pressure 0, . 

Formally, one can get phantom by switching the sign 
of kinetic energy of the standard scalar field Lagrangian, 
i.e., by taking £ = — (l/2)<9 M 0i9 M 0— 1/(0) which gives the 
energy density g = — (l/2)<j) 2 + V((f)) and the pressurep = 
-i\i^_-V^) and leads to g+p = -<j> 2 < (HHIIl 
lillll^TlMlllil20| . Phantom type of matter may also arise 
from a bulk viscous stress due to particle production pH 
or in higher-order theories of gravity |22j. Brans-Dickc 
and non- minimally coupled scalar field theories pjjj . 

The cosmological models which allow for phantom 
matter appear naturally in the mirage cosmology of 
the braneworld scenario |24| and in kinematically-driven 
quintessence (k-essence) models pB| . 

It has been shown that the matter which violates the 
strong energy condition allows the cyclic (oscillating, 
non-singular) universes 29] . It seems that similar so- 
lutions should also appear for the phantom matter which 
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violates the weak energy condition. 

A different idea of a cyclic universe in which the con- 
tracting big crunch phase can be connected with an ex- 
panding big bang phase 26] has been revived recently in 
the context of M-theory cosmology motivated ekpyrotic 
scenario [22j • In this reference it has been shown in Ref . 
|2Sj that the only way to make a transition from a con- 
tracting phase to an expanding phase in a flat universe 
is to violate the weak energy condition and this gives 
another motivation for studying phantom cosmologies. 

Some other proposals which may have something in 
common with phantom are Cardassian models in which 
one adds an extra power in the energy density term to the 
Friedmann equation [32, Om[, or modified Einstein- 
Hilbert action models in which an extra term of arbi- 
trary (both positive and negative) power of scalar cur- 
vature appears in the gravitational action j^, l3rl Is?) ]. 
Both these possibilities may be directly motivated by 
string/M-theory (HHH]. 



II. BASIC SYSTEM OF EQUATIONS WITH 
PHANTOM 



Phantom is a new type of cosmological fluid which has 
a very strong negative pressure which violates the weak 
energy condition uM > Le -' ^ obeys the equation of 
state 



P = il - 1) Q = W Q, 
with negative barotropic index 

7 = w + 1 < 0, 



(1) 



(2) 



and this implies that p < —g. 

In order to study phantom cosmologies, we start our 
discussion with the basic system of equations for isotropic 
and homogeneous Friedmann universe which reads as 



K 2 g = -A + 3— + 3— , 



n 2 p 



A 



(3) 
(4) 
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where a(t) is the scale factor, k = 0,±1 the curvature 
index, A the cosmological constant, k 2 — 8wG is the 
Einstein constant. 

Using equations ©-0 one gets the Friedmann equa- 
tion in the form 



A 



a 2 k 2 i 
a 2= Y P ~^ 2 + 3' 



(5) 



After imposing conservation law ga 37 — (3/k 2 )C 7 = 
const., one gets the Eq. © as follows 



1 (da 
~ I dt 



a 3 ^ 



A 
3"' 



(6) 



The cases which involve all types of cosmological fluids, 
i.e., with 7 = 4/3 (radiation), 7=1 (dust), 7 = 2/3 (cos- 
mic strings), 7 = 1/3 (domain walls) and 7 = (cosmo- 
logical constant) have been exactly integrated in terms of 
elliptic functions |2^. In fact, a large class of oscillating 
(non-singular) solutions have been found. 

Now we extend this discussion of the exact cosmo- 
logical solutions into the case of phantom matter with 
7 = —1/3 (we will call it phantom) and 7 = —2/3 (we 
will call it superphantom). One should emphasize that 
the "border" between "standard" and phantom models 
is given by the value of barotropic index 



7 = 



-1 



(7) 



Due to the appearance of the higher than four powers of 
the scale factor in the Friedmann equation one cannot 
integrate a general case which involves phantom matter 
by elliptic functions. However, a lot of interesting special 
solutions can be found without using the elliptic func- 
tions and this is the matter of the present paper. We 
also include stiff- fluid matter with 7 = 2. 
From the conservation law we have 



g oc a 



-3 7 



for 



7 > 



and 



g cx a 



3|7l 



for 7 < (phantom). 



(8) 



(9) 



From @ it is clear that big-bang and big-crunch sin- 
gularities appear at infinite values of the scale factor for 
phantom cosmologies. This is also obvious after studying 
the simplest solution of JfjJ for k = A = 0, which reads 



so that 



a(t) cx t 3 ~i , 



g oc t 



(10) 



(11) 



In other words, taking 7 = —1/3 (phantom) one has 
a(t) — > 00 and g — > 00 if t — > 0, while a(t) — > and g — ► 
if t — > 00. On the other hand, in a standard case 7 = 1/3, 
for example, one has a(t) — > and g — > 00 if t — > 0, while 
a(t) — » 00 and g — > if t — > 00. Similarly, one can 



show that the curvature invariants are proportional to 
the energy density and so they are divergent wherever 
the density /scale factor diverges 



R 



K 2 (3 7 - 2) q; 

.a r 



(37~2) 2 + i(3 7 + 2) 2 



g 2 . (12) 



Another interesting remark can be extracted from the 
equations ©-0 and ©-© if we admit shear anisotropy 
CTg/a 6 ((To — const.) and consider non- isotropic Bianchi 
type IX models. Namely, for 7 < 0, the shear anisotropy 
cannot dominate over the phantom matter on the ap- 
proach to a singularity when a — > 00, i.e., we have 



ga 



3|7l 



> 



for 



(13) 



and this prevents the appearance of chaotic behaviour of 
the phantom cosmologies of the Bianchi type IX [3(], El • 



III. PHANTOM MODELS DUALITY 

The system of equations ©-Q can be presented in 
the form of the nonlinear oscillator 



•■ D 2 , 

X AX 
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D(D- l)kX 



1-2/D _ 



= 0, 



(14) 



after introducing the variables [4l| 



X = a Di - w \ D(w) = ^(l + w). (15) 

For flat k — models the oscillator (mathematical pendu- 
lum) is in the lower equilibrium position provided A < 
and in the upper equilibrium position provided A > It 
is also easy to notice that the equation (|14fl preserves its 
form under the change 



D 



-D, 



7 — > —7 (w — > — (w + 2)) . 



(16) 



(17) 



It appears that there is a duality (similar to the scale- 
factor duality in pre-big-bang models which is motivated 
by superstring theory duality symmetries [42l l43j ) be- 
tween the scale factor and its inverse a — > a -1 for stan- 
dard matter (7 > 0) and phantom matter - 7 < models. 
Namely, the standard models for 1/7 = 7 — 1 are dual to 
phantom models for w' = — (2 + w) = —7 — 1 with re- 
spect to 7 = line. For example the domain wall models 
Tu, = 1/3 are dual to phantom models 7 = —1/3 and 
cosmic string models (7 CS = 2/3) are dual to superphan- 
tom models. However, this simple duality is valid only 
for flat models (see also the Refs. 0,23 where the sim- 
ilar results for both flat and non-flat models have been 
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obtained and in the non-flat case the analogy might per- 
haps be the non-abelian duality 46] ) . In the former case 
we have (for walls, phantom and A < 0) 



sin L^il I A It 



a ph = sin - 



D \i 

^7T 



A I t 



so that we have 



l ph ■ 



1/2 



-D 



ph • 



(18) 
(19) 

(20) 



From these we can conclude that standard fluids like: 
dust 7 = 1, radiation 7 = 4/3, stiff- fluid 7 = 2 are dual 
to phantoms with 7 = — 1 (p = — 2g), 7 = 4/3, 7 = — 2 
(p = — 3g), respectively. 

It is important to notice that duality in scale factor 
does not lead to the avoidance of singularity in the en- 
ergy density which results from ©-@ an d shows that 
whatever behaviour of the scale factor the density di- 
verges leading to a big-bang singularity. 



IV. PHANTOM COSMOLOGICAL MODELS 



A. Negative pressure fluids only phantom models 

General solutions of Eq. (|27|l are given in terms of ellip- 
tic or hyperelliptic functions. Here we concentrate only 
on exact elementary solutions leaving a general discus- 
sion for a separate paper |47l |. Firstly, we assume non- 
vanishing negative pressure fluids only (except negative 
A-term which gives positive pressure), i.e., fl m o = Qro = 
Qsto — 0, so that we have 



dy 
du 



ttspoy 4 + ttphoy 3 + ^aou 2 +fl w0 y+fLK>o- (28) 



Note that after a change of variables 

1 

P= ~, 

y 

we have 



(29) 



^ ) = Q sp0 + fl pho p+QAoP 2 + ^woP 3 + ^K'oP A , (30) 



which may also be useful for exact integration. 

Now we consider some special cases and refer to their 
mathematical discussion presented in the Appendix A. 



For the sake of a possible comparison with observa- 
tional data we introduce dimensionless density parame- 
ters m E3 



Qko — 



0x0, 



n 



A 



3#o 2 
K 

H Q a Q 

Ap 



(21) 
(22) 
(23) 



where x = sp,ph,w,cs,m,r, st (superphantom, phan- 
tom, domain walls, cosmic strings, dust, radiation and 
stiff-fluid, respectively), H = a/a, and q — —da/a 2 , while 

m 



Qk>0 — ^r.sfl — ^ 



(24) 



which means that the Friedmann equation © can be 
written down in the form 

^spO + QphO + ^A + &w0 + &<K'0 + ^mO + + ^stO = 1 • 

Now using the variables 



y 



a Q 



u = Hot, 



(26) 



one turns the Eq. 10 into the form 



dy 



-j^J = Qspoy + Qphoy + QaoV + Qwoy + Qk>o 
+ n^y- 1 + n r0V - 2 + n stQ y^ = Q(y). (27) 




FIG. 1: The "bounce" solution $T}. The big-bang/big- 
crunch singularities appear for y — » oo since g <x y — *• oo 
at these points so that this "bounce" in scale factor does not 
lead to singularity avoidance. 



1. No-wall models with Q, w q = fix'o = 

This case falls under a general classification given in 
the Appendix A with ai = IIao G R, a\ = fl p ho > 0, 
and ao = £l S pO > (we exclude the possibility of both 
phantoms lacking, i.e., fl sp o — Qpho = 0). The possible 
solutions are thus (cf. Appendix A): 

1. n sp o = 0, fl p ho > 1, and J7ao < (case 1.3.2 - cf. 

Fig. 0). 

2. Q sp o > 0, fl p ho > 1 — Q sp o, an d SIao < (case 1.3.3). 
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3. 0, sp o > 0, flpho £ 
^ao > (case 2.1). 



and 



4. fl sp Q < 1, and Oao > (case 2.2.3). 



5. fispo = 0, < Op/io < 1, and fl\o > (case 2.3.4). 



6. Q sp0 > 0, < Qpho < 1 - fi s 
2.3.5). 



and r^Ao > (case 



Most of these solutions (cases 3, 4, 5, 6) pass through 
zero in x, so there arise infinite values of y for finite times. 
Bearing in mind Eq. ||§J one can notice that this corre- 
sponds the energy density approaching zero. However, 
in the cases 1 and 2 one obtains cyclic "bounce" in scale 
factor y solution (which is also a "bounce" in the energy 
density g but not singularity avoidance since g — > oo at 
finite u) of the form: 

2fiA0 



^Im ~ 4fi sp o^Ao sin[i/|f2Ao|(w - uo)] - &pho 

.(31) 

which is shown in Fig. ^ In fact, there is a competition 
between the positive pressure of a negative cosmological 
term and the negative pressure of a phantom, but one is 
not able to avoid singularity on the same basis as with 
negative cosmological constant and domain walls |2flj . 

Similar behaviour is present in the case 3, though now 
there is only one infinity - either the scale factor collapses 
from infinite into a finite size in finite time, and continues 
to shrink asymptotically to zero, or it expands from zero 
at u = — oo, and reaches infinite size in a finite time. This 
is shown in Fig. [3 In fact, this behaviour is similar to 
what one has in pre-big-bang cosmology |43J with y_ as 
a pre-big-bang branch and as a post-big-bang branch. 

The cases 4, 5 and 6 are analogous to the case 3. The 
appropriate formulas are obtained from those of 2.3 in 
Appendix A, but the graphs are practically like those in 
FigE 




FIG. 2: Phantom with negative Q.ao models (case 2 of Section 
IV. A. 1) with two branches separated by curvature singularity 
(compare |43jl 1 



2. Phantom only models: Oao = !i»o = ^k'o = 

If both phantoms only are present in l|28[) we have 

1 



V 



(32) 



for Qpho + Q sp o — 1 an d Q p ho 7^ (see Fig. |3| while for 
Qpho = we have 



1 

y± = ± , 

U — Uo 

and the behaviour of the model is as in Fig. [21 



(33) 





FIG. 3: The model l|32jl for two phantoms only. Two cases 
are shown: f2 3p o = 0, £l p ho = 1; ^ls P o = Opfco = 1/2. 



B. Even powers of the polynomial Q(y) only models 

In this subsection we consider only the even powers of 
the polynomial Q(y) in (|27|l . i.e., 

— J = Qspoy 4 + ^Aoy 2 + Qk>o + ^ro2T 2 + Q s toy~ i 

(34) 

and make the substitution 



1 , 2du 
z=-s, dr]= , 

y y 



(35) 



so that reads as 
dz\ 2 

—j = n sp0 + Q A0 z + n K i Q z 2 + Q r0 z 3 + Q sto z 4: (36) 

From now on let us assume that there is no radiation 
and stiff-fluid in Q34|). i.e., fl r o — Q st o = 0. The polyno- 
mial coefficients are then (cf. Appendix A): ai = Qk'o, 
ai = f^AO) an d a — £l sp o. The possible solutions are: 

1'. FIk'q < 0, CIao > 1 — fispo (case 1.3.3). 

2'. a K , > o, n A0 g 

(case 2.1). 



5 



3'. Qk'o > 0, fUo 
1, 

and f^Ao = - 
(case 2.2.1). 



= -2- v /f2 sp0 (l + \A^po) wnen ^spo < 
-2 v /r2 sp o(V^spO ± 1) when n sp0 > 1 



4'. f^'o > 0, fi sp0 < 1 and f2 A0 = 2^f2 sp0 (l - v^po) 
(case 2.2.3). 

5' n K 'o > q_ A0 < A il A <£ 

[— 2-y/f2«po(l + y^spo)j 2-y/f2spo(l — y / n^Jo)] 
(case 2.3.1). 

6' Ox/ > 0, < fi A o < 1 - ^ S)9 o (case 2.3.5). 

All of these contain cases where z passes through zero, 
and that requires a closer analysis. Firstly, because 
the solutions themselves change — we are interested in 
y(u) — 1/y z(j](u)), which, like before, may become in- 
finite. And secondly, we are now working in a modified 
conformal time rj, and to analyse the solution in the cos- 
mological time t = u/H Q , we need consider the conver- 
gence of the integral 



Fortunately, all the solutions, are of trygonometric or 
exponential form, so that the integral reduces to a con- 
vergent elliptic one. In other words, the values of u are 
finite for finite values of r\. Moreover, in some cases the 
integral converges with r\ — ► oo. The details are given in 
particular cases. 

Although the relation between y and z is different than 
in subsection IV. A, the qualitative properties remain the 
same, that is to say, where the same classes of solutions 
apply. 

Thus, the case V is again a "bounce" shown in Fig. ^ 
The case 2' is a peculiar model in which the aforemen- 
tioned integrals converge, and u has finite values for both 
the point in which the size is infinite, and when it is zero 
(rj — > oo). This is depicted in Fig. ^ 




FIG. 4: The model with superphantom, cosmological constant 
and fVo / (solution IV.B.2'). 

For the solutions of the case 3' the shape of the func- 
tion z(rj) introduced in (|35|l is given in Fig[5] of Ap- 
pendix A with x replaced by z and r\ replaced by u in 



the graph. However, we also present these solutions in 
terms of the scale factor y(u) in Fig|S] From the dia- 
gram we can see that there is a static solution y = const, 
which corresponds to a double root of equation (|34|l with 
f2 r o = fi s to = 0. Note that the static model falls outside 
the classification given in the Appendix A since the Hub- 
ble parameter is equal to zero in this case (which is the 
result of the rescaled definition of the time parameter u in 
(|26[n . Apart from the static model we have four asymp- 
totic solutions: two of them asymptotically approach the 
static model in future infinity while the other two start 
asymptotically with the static model at past infinity. 



FIG. 5: The solution of the case B.3'. Apart from the static 
model (straight line in the middle) there are four asymptotic 
solutions: two of them asymptotically approach static model 
in future infinity while the other two start asymptotically with 
the static model at past infinity. 

The cases 4' and 6' are similar to the case 2' shown in 
Fig. H 

For the solutions of cases 5', the shape of the function 
z(rj) introduced in l|35|l is given in Fig llOl of Appendix A 
with x replaced by z and r\ replaced by u in the graph. 
In terms of the scale factor y(u) these solutions are given 
in Fig|SJ One of them is similar to that of Fig^ (a cycle 
from y cx g — > oo to a minimum in y and g and again to 
y cx q —>■ oo) and the second describes a cycle in y and g 
from zero to a maximum and again to zero. 



C. Phantom, walls and A-term models only 

This is, in fact, less general case than Case B.l, so it is 
a straightforward task to apply those solutions here. We 
have 



dy - 

3- = M pho y + tt A0 y + fl wQ y + U K > + Vt m0 y , 



which after the change of variables 

1 , du 
p=-, dr]= — 

y Vy 



(37) 
(38) 
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gives 



dy 



dr) 



ttphO + flAaP + tliuOp 2 + Qk'OP 3 + ftmOP 4 , (39) 



and this equation is integrable in terms of elliptic func- 
tions. 



D. Radiation and superphantom only models 

Leaving only the radiation pressure £l r o and super- 
phantom fl S pQ non- vanishing in l|27|) and making the sub- 
stitution 



y 



drj = y 4 du, 



we get the simple equation 



(dv \ 2 
J = Qspor + fi r o, 



which solves by 



i + n 



spO 



spO 



(40) 



(41) 



(42) 



and we have used the condition that fi r o = 1 — Q S pO- 

The typical evolution as considered in terms of the 
scale factor y is similar to that given in Fig^l Despite the 
fact that the scale factor y reaches either zero or infinity 
in both situations there are singularities - for y — > it 
is dominated by the radiation so g — > infty while for 
y — > oo it is dominated by superphantom and g — > oo, 
too. 




FIG. 6: The solutions of the case B.5'. The first is a cyclic 
solution from y oc g — » oo to a minimum in y and g and again 
to y oc g — > oo) and the second is a cycle in y and g from zero 
to a maximum and again to zero. 



E. Dust and phantom only models 

Although this case is not elementary, we present it as 
an important result in which we do not abandon the basic 



component of the universe. Takin all the other terms in 
(|27|l equal to zero except for dust f2 m0 and phantom fl p ho 
and making the change of the time coordinate u as follows 



we have the equation 



4 . ^mG 



phO 



which solves by 



4fW 2 (?7)-ft m0 ' 



(43) 



(44) 



(45) 



with V(r]) is the Weierstrass elliptic function |29| . 

Despite this solution being elliptic, the evolution is also 
as in Fig 01 The first singularity appears when y — > 
and it is dominated by dust while the second singularity 
appears when y — > oo and it is dominated by phantom. 



V. CONCLUSION 

We have studied cosmological models with phantom 
p < —g matter which violates the weak energy condition. 
We have shown that phantom matter characterised by 
negative barotropic index 7 < allows the curvature 
singularity for infinite values of the scale factor both in 
the past and in the future. We have also shown that 
the singularity cannot be dominated by shear so that 
unlike for standard 7 > cosmologies in a class of non- 
isotropic Bianchi type IX phantom models there should 
be no chaotic behaviour on the approach to a singularity. 
We have discussed a simple model of duality of the scale 
factor for phantom (7 < 0) and standard (7 > 0) matter 
solutions which is alike the scale factor duality of pre-big- 
bang cosmology motivated by superstring theory duality 
symmetries. 

We have presented a series of exact phantom cosmolo- 
gies which integrate elementary, leaving the investigation 
of the full class of models which can be in tegr ated in 
terms of elliptic functions for a future paper |42|. Apart 
from phantom matter we have admitted most of the 
known types of matter content in the universe with dis- 
crete barotropic index 7 in the equation of state such as: 
dust, radiation, stiff-fluid, domain walls, cosmic strings, 
and cosmological constant. 

One interesting class of solutions we obtained contains 
models which start with singularity in which the scale 
factor is infinite (a — > 00), then decreases, reaching some 
minimum a = a m i n , finally expanding again into another 
singularity where (a — * 00). In the standard picture for 
7 > matter, the universe follows the sequence: expan- 
sion - maximum - recollapse. On the other hand, for 
phantom 7 < models one has the opposite sequence: 
collapse - minimum - expansion. Another interesting 
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class of solutions is characterised by either an expansion 
from zero energy density to a future curvature singularity, 
or a collapse from the infinite value of the scale factor into 
a state of zero energy density. There exists also static so- 
lutions (generalized Einstein Static Universes) together 
with asymptotic solutions though asymptotic solutions 
approach static model either from singularity of the en- 
ergy density (in which the scale factor a — > oo), or from 
the state of zero energy density (in which the scale factor 
a -> 0). 

Finally, in the models which contain a mixture of phan- 
tom and standard matter (e.g. dust, radiation) there 
exist two types of singularities (with the energy density 
g — > oo): one is dominated by phantom with the scale 
factor y — > oo and another is dominated by standard 
matter with the scale factor y — > 0. 

We remark that phantom cosmologies form an inter- 
esting set of the models of the universe which do not 
contradict observational data of supernovae type la and 
that their mathematical and physical properties deserve 
further studies. 



1. a 2 < 
1.1. A < 

This makes evolution impossible, as W(x) < for all 
real values of y. Writing A = a\ + 4a ai + 4a (a — 1), 
we get its discriminant A = 16a , and the roots: a 1=F = 
— 2^/ao(y/ao ± 1). Clearly, ao < 0, would leads to A > 0, 
and ao = implies a 2 = which contradicts l|A2(l . For 
ao > 0, we would need a\ £ (&i_; &1+), but this does not 
hold together with l|A2|l . 

Such case is impossible with the assumptions made. 



1.2. A = 

As above we could have A = 0, but this would lead 
to ao = a\ = 0, and break (|a1)| . If A > the only 
possibilities, namely a\ = a^, are ruled out by the same 
condition. 
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Here we could have ao > 0, and a\ > a\ + or a\ < ai_ 
which, reasoning as above, gives in the end 01 = 1 — 
a 2 — o-q > 1 — ao- Alternatively, ao < 0, leads directly to 
ai > 1 — oq. 



APPENDIX A: GENERAL CLASSIFICATION OF 
SOLUTIONS TO THE CANONICAL EQUATION 

We take into account the equation 



^) = a 2 X 2 +.7j,r -!- (In = U'(.f) (All 



together with the constraint: 



1.3.1. ei > e 2 > 

This requires a\ > 0, ao < which can hold with to- 
gether with ai > 1 — a . The solution is periodic in the 
form: 



\/Acos[ v /ja 2 T(w - u )] - ai 

202 



(A4) 



a 2 + ai + a = 1, 



(A2) 



and seek its solutions, that might depict the evolution of 
the Universe, that is, such that x > 0. We also define e\ 
and e-i to be the roots of the quation W(x) = 0, and A = 
a\ — 4aoa2- The main constraint on y is the nonnegativity 
of the polynomial W(x) implied by equation (|A1|) . 
The general solution of the equation lAlfl is 



2a,2X + ai 



A 



= cosh[ v / a 2 "('a — u )] 



(A3) 



where uq can, in general, be complex. The graphs depict- 
ing the particular cases, were drawn with uq = 0, except 
for case 2.3.1 whose solution requires Uq = i7r/2 v /a2. It 
should be kept in mind, that the solutions may be arbi- 
trarily shifted in the direction of u-axis. 



1.3.2. ei > e 2 = 

Necessarily ao = and a\ = 1 — 0,2 > 1. The solution 
is simply: 



V = 



2a 2 



(^COs[y/\a2~\(u - U )] 



1 



(A5) 



It passes through zero, unlike the previous one. 



1.3.3. ei > > e 2 

The conditions here are: ao > and a\ > 1 — ao. The 
solution now passes through zero twice, which means it 
is no longer periodic, but depicts a single "cycle" . The 
formula is identical to (IA4I) . 
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1.3.4. > ei > e 2 

This case would require a < and a\ < 0, which 
cannot hold together with 1|A2[) . The plots of 1.3.1, 1.3.2 
and 1.3.3 are given in Fig0 



a 2 > 



When < x < ex, the expansion is asymptotical, and 
collapse occurs in finite time. 



x± = — (1 - e ± ^ u - Uo) ) 
2a 2 



(A8) 



Also, x = ei is a static solution which is unstable due to 
nearby asymptotic solutions. The plot of 2.2.1 is given 
in FigEI 



2.2.2. ei = e 2 = 




We immediately get ax 
following solution: 

x± 



ao = which yields the 



FIG. 7: Cases 1.3.1, 1.3.2 and 1.3.3 



2.2.3. ei = e 2 < 

ai = ai+ > 0, requiring ao < 1. This case is almost 
identical to 2.2.1, only the root is now negative, eliminat- 
ing the solutions x < e\. The only possibility for collapse 
is to happen in finite time now. 



2.1. A < 

This is only possible with A > => ao > 0, oi £ 
(ai_; di+), which is now a stronger restriction than a\ = 
1 — a-i — ao < 1 — ao- There are two solutions now, 
depending on the choice of the derivative sign in equation 



x± 



±y|A|sinh[ % /ai'(u - u )] - a i 



2a 2 



(A6) 



These are respectively, monotonic expansion (from zero 
size to infinity in infinite time) or monotonic collapse (in 
finite time from finite size). This is plotted in FigO 



2.2. A = 

Note, that ao < would lead to A < and hence 
A > 0, and this is not the case. 



2.2.1. ei = e 2 > 

ao > 0, so ai = a\- < 0, and for ao > 1 we also have 
ax = ax+ < 0. For x > ex = the evolution is either 
an unbounded expansion, or an asymptotical collapse. 



x± 



-ax 
2a 2 



(1 + e ± V S 2( !i ~ tl 0)-) 



(A7) 




FIG. 8: Case 2.1 




FIG. 9: Case 2.2.1. Case 2.2.2 may be pictured if a\ — 0; case 
2.2.3 when -ai/2a 2 < 0. 
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2.3. A > 



2.3.2. ei > e 2 = 



The condition a± < 1 — ao applies to all the subcases 
here. Depending on ao there are also further restrictions: 

a > => A > 0, so oi <^ [oi_;oi+]. 

ao = ai 7^ 0. 

a < =>■ No other conditions. 

In general, the situation resembles that in 1.3. All the 
graphs look the same the only difference being a shift in 
the y direction. 



This is almost identical to the previous case, only now 
x_ passes through zero, which is its only nonnegative 
value, thus giving rise to a static universe of zero size. 
x + behaves as above. Here, we have: ao = and a\ < 0. 



2.3.3. ei > > e 2 

ao < and a\ < 1 — ao now. Again, it is like 2.3.1 with 
_ discarded. 



2.3.1. ei > e 2 > 

ao > 0, and > a\ give two admissible regions of 
evolution. For x > ei, we have a "bounce" (infinite size 
is reached in infinite time, though). When < x < ei, 
the evolution is at most one "cycle" similar to that of 
1.3.3. To be exact: 



x± 



±VA cosh[ v / 02(w — uq)] — ai 

202 ' 



(A10) 



The plot of 2.3.1 is given in FiglTOl 



2.3.4 ei = > e 2 

ao = and < ai < 1. This is another "bounce" 
solution which passes through a possible singularity at 
x = 0. 



2.3.5. > ei > e 2 

ao > and < a\ < 1 — ao. The solution is sepa- 
rated into two independent ones, describing expansion or 
collapse, both going through zero. 
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